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ABSTRACT 



Using archival as well as freshly acquired data, we assess the X-ray be- 



haviour of the 



log (lAbdo et al. 



ermi/ 



jAT 7-ray pulsars listed in the First Fermi source cata- 
2010cl ). After revisiting the relationships between the pulsars' 
rotational energy losses and their X and 7-ray luminosities, we focus on the 
distance-indipendent 7 to X-ray flux ratios. When plotting our F^/Fx values as 
a function of the pulsars' rotational energy losses, one immediately sees that pul- 
sars with similar energetics have F^/Fx spanning 3 decades. Such spread, most 
probably stemming from vastly different geometrical configurations of the X and 
7-ray emitting regions, defies any straightforward interpretation of the plot. In- 
deed, while energetic pulsars do have low F^/Fx values, little can be said for 
the bulk of the Fermi neutron stars. Dividing our pulsar sample into radio-loud 
and radio-quiet subsamples, we find that, on average, radio-quiet pulsars do have 
higher values of F^/Fx, implying an intrinsec faintness of their X-ray emission 
and/or a different geometrical configuration. Moreover, despite the large spread 
mentioned above, statistical tests show a lower scatter in the radio-quiet dataset 
with respect to the radio-loud one, pointing to a somewhat more constrained 
geometry for the radio-quiet objects with respect to the radio-loud ones. 



Subject headings: gamma rays: general, x rays: general, pulsars: general, stars: 
neutron 
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Introduction 



The vast majority of tlie 1800 rotation-powered pulsars known to date (IManchester et al. 



20051 ) were discovere d by radio 



optical band (see e.g. 



Mignani 



telescopes. While only few pulsars have also been seen in the 



2008 



20101 ). the contribution of Chandra and XMM- Newton 



telescopes increased the n umber o : 



total of to ~100 (see e.g. 



Becker 



' X-r ay counterparts of radio pulsars bringing the gran 



20091 ). Such high-energy emission can yield crucial 
information on the pulsar physics, disentangling thermal components from non-thermal 
ones, and tracing the presence of pulsar wind nebulae (PWNe). 

Chandra's exceptional spatial resolution made it possible to discriminate clearly the 
PWN and the PSR contributions while XMM-Newton's high spectral resolution and 
throughput unveiled the multiple spectral components which characterize pulsars (see e.g. 



Possenti et al. 



( I2OO2I )). Altough the X-ray n on-thermal power 



Kaspi et al. 



aw index seems somehow 



2OO4J ). extrapolating the X-ray 



related to the gamma-ray spectrum (see e.g. 
data u nderpredicts the 7-ray flux by at least one order of magnitude (see e.g. 
2010d[ ). 



Abdo et al. 



Until the launch of Fermi, only seven pulsars were seen in high-energy gamma rays 



( iThompson 



20081 ) ■ and only one of them, Geminga, was not detected by radio telescopes. 
The Fermi Large Area Telescope (LAT) changed dramatically such scenario estabilishing 
radio-quiet pulsars as a major family of 7-ray emitting neutron stars. After one year 
of all-sky m onitoring Fermi/LAT has detected 54 gamma-ray pulsars, 22 of which are 



radio-quiet (Abdo et al. 



2010a 



Saz Parkinson et al 



2010; 



Camilo et al.l 1200911 . Throughout 



this pape r we shall c 



searches ( lAbdo et al. 



assify as radio-quiet all the pulsar s detected by Fermi through blind 



2010a 



Saz Parkinson et al. 



20101 ) but not seen in radio in spite of 



dedicated deep searches. Containing a sizeble fraction of radio-quiet pulsars, the Fermi 
sample provides, for the first time, the possibility to compare the phenomenology of 
radio-loud and radio-quiet neutron stars assessing their similarities and their differences (if 
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any). 

While our work rests on the Fermi data analysis and results (lAbdo et al 



Saz Parkinson et al 



2010a : 



20101) for the X-ray side we had first to build an homogeneus data set 
relying both on archival sources and on fresh observations. 

In the following we will address the relationship between the classical pulsar parameters, 
such as age and overall energetics E, and their X and 7-ray yields. While the evolution of 
the X and 7-ray luminosities as a function of E and the characteristic age have been 
already discussed, we will concentrate on the ratio between the X and 7-ray luminosities 
thus overcoming the distance conondrum which has hampered the studies discussed so 
far in the literature. We note that F^/Fx parameter probes both pulsar efficiencies at 
different wavelenghts and distribution of the emitting regions in the pulsar magnetosphere. 
Thus, such a distance-indipendent approach does magnify the role of both geometry and 
geography in determining the high-energy emission from pulsars. 



2. Data Analysis 



2.1. 7-ray Analysis 



(Abdo et al. 



We consider a ll the pulsars listed in the First Year Catalog of Fermi 7-ray sources 



(Abdo et al. 



2010a) which contains the 7-ray pulsars listed in the Firs t Fermi pulsar catalo g 



2010al ) as well as the new blind search pulsars found by 



Saz Parkinson et al 



( 120101 ). Our sample comprehends 54 pulsars: 

- 29 detected using radio ephemerides 

- 25 found through blind searches; of these 3 were later found to have also a radio emission 
and, as such, they were added to the radio emitting ones. 

Thus, our sample of 7-ray emitting neutron stars consists of 32 radio pulsars and 22 



radio-quiet pulsars. Here, we summarize the main characteristics of the analysis performed 
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in the two articles. 

The pulsar spectra were fitted with an exponential cutoff powerlaw model of the form: 
dN/dE = KE-lyexp{-E/E,utoff) 

1 GeV has been chosen to define the normalization factor because it's the energy at 
which the relative uncertainty on the differential flux is minimal. 
The spectral analysis was performed taking into account the contribution of all the 
neighboring sources (up to 17°) and the diffuse emission. Sources at more than 3° from any 
pulsars were assigned flxed spectra, taken from the all-sky analysis. 7-rays with E>100 
MeV have been used and the contamination produced by cosmic rays interactions in the 
Earth's atmosphere was avoided by selecting a zenith angle >105°. 

At flrst, all events have been used in order to obtain a phase-averaged spectrum for each 
pulsar. Next the data have been sphtted into on-pulse and off-pulse samples. The off-pulse 
sample has been described with a simple powerlaw while, for the on-pulse emission, an 
exponentially cutoff powerlaw has been used, with the off-pulse emission (scaled to the 
on-pulse phase interval) added to the model. Such an approach is adopted in order to avoid 
a possible PWN contamination to the pulsar spectrum. 

For completeness, w e included in ou r sample also the 4 radio pulsars listed in the 4th 



IBIS/ISGRI ca talog ( 



Fermi catalog (lAbdo et al. 



Bird et al. 



20091 ) but, so far, not seen by Fermi. Searching in the 1-year 



20104), we found a potential counterpart for PSR J0540-6919 



but the lack of a pulsation prevent us to associate the IBIS/ISGRI pulsar with the Fermi 
source. We therefore used the IFGL flux as an upper limit. The three remaining IBIS 
pulsars happen to be located near the galactic centre, where the intense radiation from the 
disk of our Galaxy hampers th e detection of 7-ray sources. We used the sensitivity map 



taken from I Ab do et al 



(l2010cl ) to evaluate the Fermi flux upper limit. 
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2.2. X-ray Data 

The X-ray coverage of the Fermi LAT pulsars is uneven since the majority of the newly 
discovered radio-quiet PSRs have never been the target of a deep X-ray observation, while 
for other well-known 7-ray pulsars - such as Crab, Vela and Geminga - one can rely on a 
lot of observations. To account for such an uneven coverage, we classify the X-ray spectra 
on the basis of the public X-ray data available, thus assigning: 

- label "0" to pulsars with no confirmed X-ray counterparts (or without a non-thermal 
spectral component); 

- label "1" to pulsars with a confirmed counterpart but too few photons to assess its spectral 
shape; 

- label "2" to pulsars with a confirmed counterpart for which the data quality allows for the 
analysis of both the pulsar and the nebula (if present). 

An "ad hoc" analysis was performed for seven pulsars for which the standard analysis 
couldn't be applied (e.g. owing to the very high thermal component of Vela or to the 
closeness of J1418-6058 to an AGN). Table 2 provides details on such pulsars. 
We consider an X-ray counterpart to be confirmed if: 

- X-ray pulsation has been detected; 

- X and Radio coordinates concide; 

- X-ray source position h as been validated through the blin d-search algorithm developed by 



the Fermi collaboration (jAbdo_et_al 



20091; 



Ray et al. 



2OI0I ). 



If none of these conditions apply, 7-ray pulsar is labelled as " 0" . 

According to our classification scheme we have 14 type-0, 7 type-1 and 37 type-2 pulsars. In 
total 44 7-ray neutron stars, 31 radio-lound and 13 radio-quiet have an X-ray counterpart. 

Since the X-ray observation database is continuously growing, the results available 
in literature encompass only fractions of the X-ray data now available. Moreover, they 
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have been obtained with different versions of the standard analysis softwares or using 
different techniques to account for the PWN contribution. Thus, with the exception of the 
well-known and bright X-ray pulsars, such as Crab or Vela, we re-analyzed all the X-ray 
data pubblicly available following an homogeneous procedure. If only a small fraction of 
the data are publicly available, we quoted results from a literature search. 

In order to assess the X-ray spectra of Fermi pulsars, we used ph otons with 



energy 0.3<E<10 ke V collected by Chan dra/ACIS (IGarmire et al 



(IStruder et al. 



200 If ). flTurner et al 



20031). XMM- Newton 



200 If ) and SWIFT/XRT f,Burrows et al. 



2005f ). We 



selected all the public observations (as of April 2010) that overlap the error box of Fermi 
pulsars or the Radio coordinates. 

We neglected all Chandra/HRC observations owing to the lack of energy resolution of 
the instrument. To analyze Chandra data, we used the Chandra Interactive Analysis of 
Observation software (CIAO version 4.1.2). The Chandra point spread function depends on 
the off-axis angle: we used for all the point sources an extraction area around the pulsar 
that cointains 90% of the events. For instance, for on-axis sources we selected all the 
photons inside a 2" radius circle, while we extracted photons from the inner part of PWNs 
(excluding the 2" radius circle of the point source) in order to assess the nebular spectra: 
such extended regions vary for pulsar to pulsar as a function of the nebula dimension and 
flux. 

We analyzed all the XMM- Newton data (both from PN and MOSl/2 detectors) with the 
XMM- Newton Science Analysis Software (SASvS.O). The raw observation data files (ODFs) 
were processed using standard pipeline tasks (epproc for PN, emproc for MOS data); we 
used only photons with event pattern 0-4 for the PN detector and 0-12 for the MOSl/2 
detectors. When necessay, a n accurate screening for soft proton flare events was done. 



following the prescription by 



De LucafcMolendi 



(12004) 



If, in addiction to XMM data, deep Chandra data were also available, we made an XMM 
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spectrum of the entire PSR+PWN and used the Chandra higher resolution in order to 
disentangle the two contributions. When only XMM-Newton data were available, the point 
source was analyzed by selecting all the photons inside a 20" radius circle while the whole 
PWN (with the exception of the 20" radius circle of the point source) was used in order to 
assess the nebular spectrum. 

We analyzed all the SWIFT/XRT data with HEASOFT version 6.5 selecting all the photons 
inside a 20" radius circle. If multiple data sets collected by the same instruments were 
found, spectra, response and effective area files for each dataset were added by using the 
mathpha, addarf and addrmf HEASOFT tools. 



19961 ) choosing, whenever 



All the spectra have been studied with XSPEC v. 12 ( lArnaudI 
possible, the same background regions for all the different observations of each pulsar. All 
the data were rebinned in order to have at least 25 counts per channel, as requested for the 
validity of statistic. 



The XMM-Chandra cross calibration studies ( IStuhlinger et al. 



20081) report only minor 



changes in flux (<10%) between the two instruments. When both XMM and Chandra data 
were available, a constant has been introduced to account for such uncertainty. Conversely, 
when the data were collected only by one instrument, a systematic error was introduced. 
All the PSRs and PWNs have been fitted with absorbed powerlaws; when statistically 
needed, a blackbody component has been added to the pulsar spectrum. Since PWNs 
typically show a powerlaw spectrum with a photon index which steepens moderately as 



a function of the distance from the PSR (IGaenslerfcSlane 



20061 ). we used only the inner 



part of each PWN. Absorption along the line of sight has been obtained through the fitting 
procedure but for the cases with very low statistic for which we used informations derived 
from observations taken in different bands. 
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2.3. X-ray Analysis 

For pulsars with a good X-ray coverage we carried out the following steps. 
If only XMM-Newton pubhc observations were available, we tried to take into account the 
PWN contribution. First we searched the literature for any evidence of the presence of a 
PWN and, if nothing was found, we analyzed the data to search for extended emission. 
If no evidence for the presence of a PWN was found, we used PN and MOSl/2 data in 
a simultaneous spectral fit. On the other hand, if a PWN was present, its contribution 
was evaluated on a case by case basis. If the statistic was good enough, we studied 
simultanously the inner region, containing both PSR and PWN, and the extended source 
region surrounding it. The inner region data were described by two absorbed (PWN and 
PSR) powerlaws, while the outer one by a single (PWN) powerlaw. The N^ and the PWN 
photon index values were the same in the two (inner and outer) datasets. 
When pubhc Chandra data were available, we evaluated separately PSR and PWN (if any) 
in a similar way. 

If both Chandra and XMM pubhc data were available, we exploited Chandra space 
resolution to evaluate the PWN contribution by: 

- obtaining two different spectra of the inner region (a), encompassing both PSR and PWN 
and of the outer region (b) encompassing only the PWN; 

- extracting a total XMM spectrum (c) containing both PSR and PWN: this is the only 
way to take into account the XMM's larger PSF; 

- fitting simultaneously a,b,c with two absorbed powerlaws and eventually (if statistically 
significant) an absorbed blackbody, using the same Nh] a constant moltiplicative was also 
introduced in order to account for a possible discrepancy between Chandra and XMM 
calibrations; 

- forcing to zero the normahzation(s) of the PSR model(s) in the Chandra outer region 
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and freeing the other normahzations in the Chandra datasets; fixing the XMM PSR 
normahzation(s) at the inner Chandra dataset one and the XMM PWN normahzation at 
the inner+outer normahzations of the Chandra PWN. 

Only for few weh-known pulsars, or pulsars for which the dataset is not yet entirely public, 
we used results taken from the literature (see Table [21). Where necessary, we used XSPEC 
in order to obtain the flux in the 0.3-10 keV energy range and to evaluate the unabsorbed flux. 



For pulsars with a confirmed counterpart but too few photons to discriminate the 
spectral shape, we evaluated an hypothetical unabsorbed flux by assuming that a single 
powerlaw spectrum with a photon index of 2 to describe PSR+PWN. We also assumed 
that the PWN and PSR thermal contributions are 30% of the entire source flux (a sort of 
mean value of all the considered type 2 pulsars). To evaluate the absorbing column, we 
need a distance value which can come either from the ra dio dispersion or - for radio- quiet 



pulsars - from the following pseudo-distance reported in 
d = 0.51Eu^'/F^% kpc 



Saz Parkinson et al. 



fcoioh : 



where E = E34 x lO^^ erq/g and F, 



is assumed to be 1 (IWatters et al. 



■F^ in X 10 erg /cm s and the beam correction factor 



20091 ) for all pulsars. 



Then, the HEASARC WebTools (| http://h easar c . gsf c . nasa. gov /do cs/tools . html ; ) was used 
to find the galactic column density (^h) in the direction of the pulsar; with the distance 
information, we could rescale the colun in density value of t he pulsar. We found the source 



count rate by using the XIMAGE task ( iGiommi et al. 



19921 ). Then, we used the WebPimms 



tool inside the WebTools package to evaluate the source unabsorbed flux. Such a value 
has then to be corrected to account for the PWN and PSR thermal contributions. We are 
aware that each pulsar can have a different photon index, as well as thermal and PWN 
contributions so that we used these mean values only as a first approximation. All the 
low-quality pulsars (type 1) will be treated separately and all the considerations in this 
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paper will be based only on high-quality objects (type 2). 

For pulsars without a confirmed counterpart we evaluated the X-ray unabsorbed flux 
upper limit assuming a single powerlaw spectrum with a photon index of 2 to describe 
PSR+PWN and using a signal to noise of 3. 

The column density has been evaluated as above. Under the previous hypotheses, we used 
the signal-to-noise definition in order to compute the upper limit to the absorbed flux of 
the X-ray counterpart. Next we used XSPEC to find the unabsorbed upper limit flux. 

On the basis of our X-ray analysis we define a subsample of Fermi 7-ray pulsars for 
which we have, at once, reliable X-ray data (type 2 pulsars) and satisfactory distance 
estimates such parallax, radio dispersion measurement, column density estimate, SNR 
association. Such a subsample contains 24 radio emitting neutron stars and 5 radio quiet 
ones. The low number of radio quiet is to be ascribed to lack of high quality X-ray data. 
Only one of the IBIS pulsars has a clear distance estimate. Moreover, we have 4 additional 
radio-quiet pulsars with reliable X-ray data but without a satisfactory distance estimate. 



In Tables [2ll3] we reported the gamma-ray and X-ray parameters of the 54 Fermi first 
year pulsars. We also included th e four hard X-ray pulsars taken from the "4th IBIS/ISGRI 



soft gamma-ray survey catalog" (iBird et al.ll2009l ). We use E = Att'^IP/P, tq = P/2P and 
Bi^ = 3.3 X 10^9(PP)^/2 X {10km)/{Rl), where Ri^ = cP/2n, P is the pulsar spin period, P 
its deriv ative and the stand ard value for moment of inertia of the ne utron star 1= 10^^ q/c m^ 



see e.g. 



Steiner et al 



Saz Parkinson et al 



2OIOI) . Using the P and P values taken from 



Abdoetal 



fl2010ah : 



(I2OIOI ). w e computed t 



distance values are taken from 



Abdo et al. 



le values reported in Table 



(j2010a|) 



11 Most of the 



Saz Parkinson et al. 



(120101 ) (see Table 
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3. Discussion 



3.1. Study of the X-ray luminosity 



The X-ray luminosity, Lx, i s correlated with the pulsar spin-down luminosity E 



The scaling was firstly noted by 



SewardfcWand (119881 ) who used Einstein data of 22 



pulsar - most of t 
and logE. Later, 



lem just upper liraits - to derive a linear relation between \ogFf^2-4k 



eV 



BeckerfcTrumpen (119971 ) investigated a sample of 27 pulsars by using 



ROSAT, yielding the simple scaling L^^ 2.4fcey ^ rjj^g uncertainty due to soft X-ray 

absorption translates into very high flux errors; moreover it was very hard to discriminate 



between the t 



herma 



Possenti et al. 



and powerlaw spectral components. A re-analysis was performed by 
(I2OO2I ). who studied in the 2-10 keV band a sample of 39 pulsars observed 
by several X-ray telescopes. However, they could not separate the PWN from the pulsar 
contribution. Moreover, they conservatively adopted, for most of the pulsars, an uncertainty 
of 40% on the distance values . A be tter comparison with our data can be done with the 



results by iKargaltsevfcPavlovl (j2008l ). who recently used high-resolution Chandra data in 
order to disentangle the PWN and pulsar fluxes. Focussing just on Chandra data, and 
rejecting XMM observations, they obtain a poor spectral characterization which translates 
in high errors on fluxes. They also adopted an uncertainty of 40% on the distance values 
for most pulsars. Despite the big uncertainties, mainly due to poor distance estimates, all 
these datasets show that the Lx versus E relation is quite scattered. The high values of the 
Xred seeni to exclude a simple statistical effect. 



We are no w facing a different 



has improved (lAbdo et al. 



2010a 



panorama, since our abilit y to evaluate pulsars' distances 



Saz Parkinson et al 



2OIOI ) and we are now much better 



in discriminating pulsar emission from its nebula. The use of XMM data makes it possible 
to build good quality spectra allowing to disentangle the non-thermal from the thermal 
contribution, when present. In particular, we can study the newly discovered radio-quiet 
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pulsar population and compare them with the "classical" radio- loud pulsars. We investigate 
the relations between the X and 7 luminosities and pulsar parameters, making use of the 
data collected in Tables [T]|2}13i 



Using the 29 Fermi type 2 pulsars with a clear distance estimate and with a 
well- constrained X-ray spectrum, the weighted least square fit yields: 



logioL^^ = (l.ll^[!i^) + (1.04 ± 0.09)/o(7io^34 



where E = x 10^ erg/ s and Lx = L^g x lO^erg/s. All the uncertains are at 90% 
confidence level. We can evaluate the goodness of this fit using the reduced chisquare value 
xled — 3-7; a double linear fit does not significantly change the value of Xred- more precise 
way to evaluate the dispersion of the dataset around the fitted curve is the parameter: 

where ylg^ is the actual i*'* value of the dataset (in our case logioL-^^) and the expected 
one. A lesser spread in the dataset translate into a lower value of W"^. We obtain 
W"^ = 0.436 for the — E relationship. Such high values of both and Xred 
indication of an important scatter ing of the Lx values a round the fitted rela t ion. 



Our results are in agreement with 



Possenti et al. 



fl2002h : iKargaltsevfcPavlovl moj ). 



3.2. Study of the 7-ray luminosity 



The gamma-ray luminosity, L^, is correlated with the pulsar 
E. Such a trend is expe cted in many theorical models (see e.g. 



spin-down luminos ity 



MuslimovfcHardine 



pulsars (lAbdo et al. 



Zhang et al. 



2004 



2003r) a nd it's shortly discussed in the Fermi LAT catalog of gamma-ray 



2mm 



Selecting the same subsample of Fermi pulsar used in the previous chapter to assess 
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the relation between and E, we found that a hnear fit: 



log.oLl^ = (0.451°:??) + (0.88 ± 0.07)/o(7io^34 



(2) 



yields an high value of xled = 4.2. 



Inspection of the distribution of residuals lead us to try a double-linear relationship, which 
yields: 



log.oLj^ = (2.45 ± 0.76) + (0.20l°JJ)/o(7io^34 , E > E,„t 
logioLl^ = (0.52 ± 0.18) + (1.43l°:i)/o5io44 , E < E,rit 



(3a) 
(3b) 



with Er 



S.72tl'AA X W^^erg/s and Xred = 2-2. An f-test shows that the probability for 



^crit — <->• ' ^—3.44 

a chance improvem ent is 0.00011. Such a result is in agreement with the data reported 



m 



Abdo et al.l (12010aj ) for the entire dataset of Fermi 7-ray pulsars. Indeed, the xled 



obtained for the double linear fit is better than that obtained for the Lx-E relationship. 
We obtain = 0.344 for the double linear — E relationship. Both the Xred ^^"^ 
in agreement with a little higher scatter in the Lx — E graph. A difference between the 
X-ray and 7-ray emission geometries - that translates in different values of f^ and fx - could 
explain such a behaviour. 

The existence of an Ecru has been posited from the theoretical point for different pulsar 
emission models. Revisiting the outer-gap mod el for pulsars with r < 10'' yrs and assuming 



Zhang et al 



(I2OO4I ) found a sharp boundary. 



initial conditions as well as pulsars' birth rates, 
due to the saturation of the gap size, for = E. They obtain the following distribution of 
pulsars' 7-ray luminosities: 



logioL^ = logioE + const. , E < Ecrit 
logioL^ ~ O.SOlogioE + const. , E > Ecrit 



(4a) 
(4b) 



By assuming the fractional gap size from 



ZhangfcCheneJ fll997h . they obtain E, 



■'crit 



1.5 X 10'^^ P^^^er g / s. While Equation 4 is similar to our double linear fit (Equation 3), the 
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Ecrit they obtain seems to be lower than our best fit value. 



20031), the break occurs at 



On the other hand, in slot-gap models (iMuslimovfcHardingl 
about lO^^erg/s, when the gap is limited by screening of the acceleration field by pairs. 
We can see from Figure [2] that radio-quiet pulsars have higher luminosities than the 
radio-loud ones, for similar values of E. As in the Lx — E fit, we can't however discriminate 
between the two population due to the big errors stemming from distance estimate. 



3.3. Study of the 7-to-X ray luminosity ratio 

At variance with the X-ray and gamma-ray luminosities, the ratio between the X-ray 
and gamma-ray luminosities is indipendent from pulsars' distances. This makes it possible 
to significatively reduce the error bars leading to more precise indications on the pulsars' 
emission mechanisms. 

Figure [3] reports the histogram of the F^/Fx values using only type 2 (high quality 
X-ray data) pulsars. The radio-loud pulsars have < F^/Fx >~ 800 while the radio-quiet 
population has < F^/Fx >~ 4800. Applying the Kolmogoroff-Smirnov test to type 2 
pulsars' F^jFx values we obtained that the chance for the two datasets belong to the same 
population is 0.0016. By using all the pulsars with a confirmed X-ray counterpart (i.e. 
including also type 1 objects) this probability increase to 0.00757. We can conclude, with 
a 3(7 confidence level, that the radio-quiet and radio-loud datasets we used are somewhat 
different. 



3.3.1. A distance indipendent spread in F^/Fx 

Figure S] shows F^/Fx as a function of E for our entire sample of 7-ray emitting NSs 
while in Fig [5] only the pulsar with "high quality" X-ray data have been selected. Even 
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neglecting the upper and lower limits (shown as triangles) as well as the low quality points 
(see Figure E]), one immediatly notes the scatter on the F^/Fx parameter values for a 
given value of E. Such an apparent spread cannot obviously be ascribed to a low statistic. 
An inspection of Figure H] makes it clear that a linear fit cannot satisfactory describe the 
data. In a sense, this finding should not come as a surprise since Figure H] is a combination 
of Figures [T] and [2] and we have seen that figure [2] requires a double linear fit. However, 
combining the results of our previous fits (Equations 1 and 3) we obtain the dashed line in 
Figure m clearly a very poor description of the data. For E ~< 5 x 10^^ the F^/Fx values 
scatter around a mean value of ~1000 with a spread of a factor about 100. For higher E 
the values of F^/Fx seem to decrease drastically to an average value of ~50, reaching the 
Crab with F^/Fx~ 0.1. 

The spread in the F^/Fx values for pulsars with similar E is obviously unrelated to 
distance uncertainties. Such a scatter can be due to geometrical effects. For both X-ray 
and 7-ray energy bands: 

L^,x = ^nf^^xFobsD^ (5) 

where fx and f^ account for the X and 7 beaming geometries (which may or may not be 
related). If the pulse profile observed along the line-of-sight at ( (where (e is the Earth 
line-of-sight) for a pulsar with magnetic inclination a is F(a,C,0), where is the pulse 
phase, than we can write: 

2 J F{aXE, 



where / depends only from the viewing angle and the magnetic inclination of the pulsar. 
With an high value of this correction coefficient, the emission is disfavoured. Obviously 
F^/Fx=L^/Lxxfx/f7. Different f-y/ fx values for different pulsars can explain the scattering 
seen in the F^ /F x-E relationship. 



Watters et al 



(120091 ) assume a nearly uniform emission efficiency while 



Zhang et al 



(I2OOJ) 
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compute a significant variation in the emission efficiency as a function of tlie geometry of 
pulsars. In botfi cases, geometry plays an important role through magnetic field inclination 
as well as through viewing angle. 

The very important scatter found for F^/Fx values is obviously due to the different 
geometrical configurations which determine the emission at different wavelength of each 
pulsar. While geometry is clearly playing an equally important role in determining pulsar 
luminosities, the F^/Fx plot makes its effect easier to appreciate. 

The dashed line in Figures H] and E] is the combination of the best fits of L^-E and 
Lx-E relationship, considering f^=l and fx=l so that represent the hypothetical value of 
Fj/Fx that each pulsar would have if f.y=fx'- all the pulsars with a value of F^/Fx below 
the line have fx <f'y We have seen in Section [373] that the radio-quiet dataset shows an 
higher mean value of F^/Fx- This is clearly visible in Figure |5] where all the radio-quiet 
points are above the expected values (dashed line) so that all the radio-quiet pulsars should 
have fx >f'y Moreover, the radio-quiet dataset shows a lower scatter with respect to the 
radio- loud one pointing to more uniform values of f^/fx for the radio-quiet pulsars. A 
similar viewing angle or a similar magnetic inclination for all the radio-quiet pulsars could 
explain such a behaviour (see Equation 6). 

Figure E] shows the F^/Fx behaviour as a function of the characteristic pulsar age. 
In view of the uncertainty of this parameter, we have also built a similar plot using 
"real" pulsar age, as derived from the associated supernova remnants (see Figure 6). 
Similarly to the E relationship, for r < 10^ years, F^/Fx values increase with age (both 
the characteristic and real ones), while for t > lO^years the behaviour becomes more complex. 
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3.4. Study of the selection effects 



There are two main selections we have done in order to obtain our sample of pulsars 
with both good 7 and X-ray spectra (type 2). First, the two populations of radio-quiet 
and radio-loud pulsars are unveiled with different techniques: using the same dataset, 
pulsars with known rotational ephemerides have a detection threshold lower t han pulsars 



found through blind period searches. In the First Fermi LAT pulsar catalog (lAbdo et al. 



2010aJ) the faintest gamma-ray-selected pulsar has a flux ~ 3x higher than the faintest 
radio-selected one. Second, we chose only pulsars with a good X-ray coverage. Such a 
coverage depends on many factors (including the policy of X-ray observatories) that cannot 
be modeled. 

Our aim is to understand if these two selections influenced in different ways the two 
populations of pulsars we are studying: if this was the case, the results obtained would 
have be en distorced. The 7 -ray selection is discussed at length in the Fermi LAT pulsar 



catalog (iAbdo et al. 



2010a!). Since the radio-quiet population has obviously a detection 
threshold higher than the radio-loud, we could avoid such bias by selecting all the pulsars 
with a flux higher than the radio-quiet detection threshold (6xl0~^ph/cm^s). Only five 
radio-loud type 2 pulsars are excluded (J0437-4715, J0613+1036, J0751+1807, J2043+2740 
and J2124-3358) with Y^/Yx values ranging from 87 to 1464. We performed our analysis 
on such a reduced sample and the results don't change significatively. 
We can, therefore, exclude the presence of an important bias due to the 7-ray selection on 
type 2 pulsars. 

In order to roughl y evaluate the selection affecting the X-ray observations, we used the 



method developed by fiSchmidt 



19681 ) to compare the c urrent radio- 



spatial distributions, following the method also used in 



AbdoetaL 



quiet a nd -loud samples' 



(l2010al ). For each object 



with an available distance estimate, we computed the maximum distance still allowing 
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detection from Dr. 



De.fit. (F-y /Fmin 



and Frrtin are taken from 



Abdo et al 



where D^.'^t comes from Table 1, the photon flux 



fl2010a[ ): 



Saz Parkinson et al 



teoioh . We limited 



to 15 kpc, and compared V, the volume enclosed within the estimated source distance, to 
that enclosed within t he maximum distance, Knax, for a galactic disk with radius 10 kpc and 



thickness 1 kpc (as in 



Abdo et al 



(|2010aD). The inferred values of < V/Vmax > are 0.462, 
0.424, 0.443 and 0.516 for the entire gamma-ray pulsars' dataset, the radio-quiet pulsars, 
millisecond pulsars and the radio-loud pulsars. These are quite close to the expected value 
of 0.5 even if < V/Vmax >^'^ is lower than < V/Vmax If we use only type 2 pulsars we 
obtain 0.395, 0.335, 0.462 and 0.419. These lower values of < V/Vmax > indicate that we 
have a good X-ray coverage only for close-by - or very bright - pulsars, not a surprising 
result. By using the X-ray-counterpart dataset, both the radio- loud and -quiet < V/Vnax > 
values appear lower of about 0.1: this seems to indicate that we used the same selection 
criteria for the two population and we minimized the selection effects in the histogram of 
Figure |3l 

We can conclude that the 7-ray selection introduced no changes in the two populations, 
while the X-ray selection excluded objects both faint and/or far away; any distortion, if 
present, is not overwhelming. 

Only if deep future X-ray observations centered on radio-quiet Fermi pulsars will fail to 
unveil lower values of F^/F^, it will be possible to be sure that radio-quiet pulsars have a 
different geometry (or a different emission mechanism) than radio-loud ones. 



4. Conclusions 

The discovery of a number of radio-quiet pulsars comparable to that of radio-loud ones 
together with the study of their X-ray counterparts made it possible, for the first time, to 
address their behaviour using a distance independent parameters such as the ratio of their 
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fluxes at X and gamma ray wavelengths. 

First, we reproduced the well known relationship between the neutron stars luminosities and 
their rotational energy losses. Next, selecting only the Fermi pulsars with good X-ray data, 
we computed the ratio between the gamma and X-ray fluxes and studied its dependence on 
the overall rotational energy loss as well as on the neutron star age. 
Much to our surprise, the distance independent F^/Fx vahies computed for pulsars of 
similar age and energetic differ by up to 3 orders of magnitude, pointing to important 
(yet poorly understood) differences both in position and height of the regions emitting 
at X and 7-ray wavelengths within the pulsars magnetospheres. Selection effects cannot 
account for the spread in the F^/Fx relationship and any further distortion, if present, is 
not overwhelming. 

In spite of the highly scattered values, a decreasing trend is seen when considering young 
and energetic pulsars. Moreover, radio quiet pulsars are characterized by higher values of 
the Fj/Fx parameter (< F^/Fx >ri^ 800 and < F^/Fx >rq^ 4800) so that a KS test 
points to a chance of 0.0016 for them to belong to the same population as the radio loud 
ones. While it would be hard to believe that radio loud and radio quiet pulsars belong to 
two different neutron star populations, the KS test probably points to different geometrical 
configurations (possibly coupled with viewing angles) that characterize radio loud and radio 
quiet pulsars. Indeed the radio-quiet population we analyzed is less scattered than the 
radio-loud one, pointing to a more uniform viewing or magnetic geometry of radio-quiet 
pulsars. 

Our work is just a starting point, based on the first harvest of gamma-ray pulsars. The 
observational panorama will quickly evolve. The gamma-ray pulsar list will certainly grow 
and this will trigger more X-ray observations, improving both in quantity and in quality 
the database of the neutron stars detected in X and 7-rays to be used to compute our 
multiwavelength, distance independent parameter. However, to fully exploit the information 
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packed in the F^/Fx a complete 3D modelling of pulsar magnetosphere is needed to account 
for the different locations and heights of the emitting regions at work at different energies. 
Such modelling could provide the clue to account for the spread we have observed for the 
ratios between 7 and X-ray fluxes as well as for the systematically higher values measured 
for radio-quiet pulsars. 

XMM data analysis is supported by contracts ASI-INAF 1/088/06/0 and NASA NIPR 
NNG10PL01I30. Chandra data analysis is supported by INAF-ASI contract n.I/023/05/0. 
It is a pleasure to thank Niel Gehrels for granting SWIFT observations of the newly 
disocvered Fermi pulsar. We also thank Pablo Saz Parkinson and Andrea Belfiore for the 
collaboration in searching for X-ray counterparts of radio-quiet pulsars. 
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Table 1:: 

a : P, P and most of the values of the distance are taken from 



Abdo et al. 



b : Ag e derived from the associ a ted SNR. Respectively taken from 



2002) 



2ooa) 



Gorenstein et al 



1 1974) 



BietenholzfcBartell COOSi ') 



Winkler et al 



These distances are taken from 



Mieliazzo et al.l 1I2OO; 



1 201031) 



Slane et al 



BockfcGvaramadzd 1I2OO: 



Saz Parkinson et al 



Kothes et al 



Saz Parkinson et al 



I gOOfK 



1 2009 ) 



Gotthelf et al 



Hales et al 



12003), 



1 20071) 



1 2010) 



Rudie et al 



RobertsfcBrogan 



i200i 



I 200a) 



|Hwang_e^_aL 



TamfcRobertd pOOSi ) 



(2001 



Thorsett et al 



Brogan et al 



2OIOI ) and are obtained under the absumption of a beam correction factor f-y = 1 for the gamma-ray emission 



cone of all pulsars. In this wray one obtains: 

d = 0.51^34 '''V^.y^ao kpc 

where E = Es4 X lO^^e^g/g and Fry = F-y, 10 X 10^ ^ '^ erg/ cm? s. See als o 
d : Respectively taken from 
: g = radio-quiet pulsars 



Campana et al 



1 20081 ) 



Kaspi et al 



1 199a ) 



Saz Parkinson et al 



Kaspi et al 



I 2OOII ) 



I 20ld) 



GotthelffcHalpernI 120091) 



Halpern et al 



1 20071 ) 



radio-loud pulsars ; m = millisecond pulsars ; i = pulsars detected by INTEGRAL/IBIS but not yet by Fermi (see 



Bird et al 



l200g|)). 

f : Only bright PWNs have been considered (with F^™" > 1/5 F^""^). The presence or the absence of a bright PWN has been valued by re-analyzing the X-ray data 
(except for the X-ray analyses taken from literature, see the following table). 
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1.82±0.03 






5.92x10 


-5 


J0218-|~4232 


2 


Lf 




4.87t?-J^ 


7.6±4.3 


1.19±0.12 






2.05x10 


-3 


J0248+6021 





S 


<9.00 


<9.00 


80= 


2 










J0357+32 


2 


C 


64+° °^ 




8.0±4.0 


2.53±0.25 











PSR Name 


X« 


Inst*" 


T?nt 

{W'^^erg/crri^s) 


(10~^^er^/cm^, 


J0437-4715 


2 


X/C 




14.30J.7 


J0534+2200 


2 


L/ 


44300±1000 


44300±1000 


J0540-6919 


2 


Lf 


568±6 


568±6 


J0613-0200 


2* 


X 


U.ZZi_Q ;^5g 


^•^^1-0.158 


J0631+1036 





X 


<0.225 


<0.225 


J0633+0632 


1 


s 


1.53±0.51 


1.53±0.51 


J0633+1746 


2 




4 07+0-09 




J0659+1414 


2 


Lf 


4 nfi+°°3 

4.UD_o.59 




J0742-2822 





X 


<0.225 


<0.225 


J0751+1807 


2 




n 44+0- 18 


n 44+0 - 18 


J0835-4510 


2* 


L/ 


65.1±15.7 


281±67 


J1023-5746 


2* 


C 


1.61±0.27 


1.61±0.27 


J1028-5819 


1 


s 


1.5±0.5 


1.5±0.5 


J1044-5737 





s 


<3.93 


<3.93 


J1048-5832 


2* 


c+x 




0.501°:?^ 


J1057-5226 


2 


c+x 


1 CI +0.02 
J--Oi_0.13 


24.5t°j 


J1124-5916 


2 


c 


9.78tl:J8 


10.90tJi2 


J1413-6205 





s 


<4.9 


<4.9 


J1418-6058 


2 


c+x 


0.353±0.154 


0.353±0.154 


J1420-6048 


2* 


X 


1.6±0.7 


1.6±0.7 


J1429-5911 





s 


<16.9 


<16.9 


J1459-60 





s 


<3.93 


<3.93 



(1020cm-3) 
4.4 ± 1.7 

34.5±0.2 
37±1 

r 

20" 
1.07'= 
4.3±0.2 
20" 

46 

2.2±0.5 

50^ 
50^ 



90 



+40 
-20 



2.7±0.2 



30.0 



+2.8 
-4.8 



40^ 
900 + 134 

9r,<j + 161 

-^"-''^-106 

80° 
100^= 



Jx 

3.17 ± 0.13 

1.63±0.09 
1.98±0.02 
2.7±0.4 

2 

2 

1.7±0.1 
2.1±0.3 
2 

1.59±0.30 
2.7±0.6 

1 ir+0.24 
l-l"J-0.22 

2 
2 

2.4±0.5 
1.7±0.1 



1.54 



+0.09 
-0.17 



i-O'J-o.se 



0.84: 



+0.55 
-0.37 



kT 
(keV) 



(km) 



0.228lg:gg| 0.0601U: 



008 



0.190±0.030 0.04±0.01 
0.125±0.003 1.80±0.15 



0.129±0.007 



2.5±0.3 



0.179±0.006 0.46±0.06 

42fi+° °3'* n O7A+0-0S9 

U.4ZD_Qoi8 ^•^'*-0.077 



Eff^ 



8.23x10" 
3.67x10- 
1.13x10- 
3.74x10- 



8.99x10- 
8.46x10- 

2.52x10- 
9.78x10- 



1.74x10- 
2.49x10- 
2.27x10- 

1.05x10- 
1.11x10- 



PSR Name 


X« 


Inst*" 


T?nt 

{W'^^erg/crri^s) 


'ptot 

{1Q~^^ erg /cm'^. 


J1509-5850 


2 


C+X 




0.89lt°;}|^ 


J1614-2230 





c+x 


<0.286 




J1617-5055 


2 




64.2+0.3 


64.2+0.3 


J1709-4429 


2 


C+X 


q 7Q + 0.37 
"J- ' *-0.94 




J1718-3825 


2 


X 


2.80+0.67 


2.80+0.67 


J1732-31 





s 


<2.42 


<2.42 


J1741-2054ff 


1 


s 






J1744-1134 





c 


<0.272 


0.272+0.020 


J1747-2958 


2* 


c+x 


48.7121^3 


48.7^21^3 


J1809-2332 


2 


c+x 


1 40+0-25 

J--*U_o.23 


o 14+0.57 
'^•J-^-0.53 


J1811-1926 


2 


c 


26.6l3:? 


26.613.7 


J1813-1246 


1 


s 


9.675+3.225 


9.675+3.225 


J1813-1749 


2 


c 


24.4+11.5 


24.4+11.5 


J1826-1256 


2 


c 


1.18+0.58 


1.18+0.58 


J1833-1034 


2 


x+c 


66.3+2.0 


66.3+2.0 


J1836+5925 


2 


x+c 


0.459t°;f,l 




J1846+0919 





s 


<2.92 


<2.92 


J1907+06 


1 


c 


3.93+1.45 


3.93+1.45 


J1952+3252 


2 




35.0+4.4 


38.0+3.0 


J1954+2836 





s 


<3.65 


<3.65 


J1957+5036 





s 


<2.98 


<2.98 


J1958+2841 


1 


s 


1.57+0.53 


1.57+0.53 



(1020cm-3) 
80^ 

345+14 
45.6t^;^ 

70*= 

50^= 

0*= 

10+42 
■■-^-12 

2561^ 



61 



+9 



175 



+11 
-12 



100^ 

100^ 
230^ 

-0.792 







-0 

20^ 



398tf7| 
30+1 



Jx 

1.31+0.15 

2 

1.14+0.06 
1.88+0.21 
1.4+0.2 
2 



2 in+°-50 

^•^^"-'-0.28 



1 "^1+012 
■•■•"3 ■'■-0.44 



1.85 



+ 1.89 
-0.36 



0.91 



+0.09 
-0.08 



2 

1.3+0.3 
0.63toj° 
1.51+0.07 



1.56 



+0.51 
-0.73 



3.16+2.28 

i.DO_Q Q5 



kT 
(keV) 



(km) 



0.236+0.024 0.92l0;^| 



r, 979+0.094 1 09 + I.6OO 



0.056 



+0.012 
-0.009 



Eff^ 



1.12x10" 



2.03x10" 



0.166+0.012 4.3lo:^6 6.62x10" 

3.12x10" 



9.93x10" 



7.41x10" 



0.190+0.025 1.54+1:11 8.98x10" 

1.18x10" 
1.13x10" 



0.13+0.02 



4.15x10" 
4.47lf:^f 5.84x10" 



2.2tli 3.57x10" 



PSR Name 


X° 


Inst'' 




T7<tot 

(lO'^^erg/cTO^s) 


(lO^Ocrn-^) 


ix 


kT 
(keV) 


(km) 


Effi 




J2021+3651 


2 


c+x 


9 r,i +0.35 




65.5+6.0 


2+0.5 


0.140tS;Sf8 


4.94+1.40 


2.75x10" 


-5 


J2021+4026 


1 


c 


0.443+0.148 


0.443+0.148 


40^= 


2 


- 


- 


1.03x10- 


-4 


J2032+4127 


2* 


c+x 


0.423+0.118 


0.423+0.118 


qo 7+75.6 
"3°- '-38.7 


1 07+0-96 


- 


- 


1.99x10- 


-4 


J2043+2740 


2 


X 


0.208l^i«° 


0.208t'i';o^ 


r, + 20 


3.1+0.4 


- 


- 


1.44x10- 


-4 


J2055+25 


2 


X 


0.3821°:!^^ 


0.3821^:1^^ 


7 q+10.4 
' -"^-7.3 


9 9+0.5 
^■^-0.6 






1.79x10- 


-3 


J2124-3358 


2 


X 




0.9591°^!^ 


2.76l^;«^ 


2 89+"-'^^ 

^•°''-0.35 


0.268i°:g^^ 


u.uia^o.oog 


1.25x10- 


-4 


J2229+6114 


2 


c+x 






3011 


1 01+"°^ 

^•'J^-0.12 






2.90x10- 


-4 


J2238+59 





s 


<4.49 


<4.49 


70^ 


2 











Table 2:: X-ray spectra of the pulsars. The fluxes are unabsorbed and here the non-thermal and total fluxes are shown. The model used is an absorbed powerlaw 
plus blackbody, where statistically necessary. The only exceptions are PSR J0437-4715 (double PC plus powerlaw), J0633+1746 and J0659+1414 (double BB plus 
powerlaw): here only the most relevant thermal component is reported. All the errors are at a 90% confidence level. 

a : This parameter shows the confidency of the X-ray spectrum of each pulsar, based on the available X-ray data. An asterisk mark the pulsars for which ad ad- hoc 
analysis was necessary. See section 12.21 

b : C = Chandra/ACIS ; X = XMM/PN+MOS ; S = SWIFT/XRT ; L = literature. Only public data have been used (at December 2010). 

c : here, the column density has been fixed by using the galactic value in the pulsar direction obtained by Webtools ( |http://heasarc.gsfc.nasa.gov/docs/tools.html^ 
and scaling it for the distance (see Table [TJ. 

d : the beam correction factor fx is assumed to be 1, which can result in an efficiency > 1. See 
e : The statistic is very low so that it was necessary to freeze the column density parameter; the values have been evaluated by using WebTools 
l |http://heasarc.gsfc.nasa.gov/docs/tools.html 



Watters et al 



I 2009 ). Here the errors are not reported. 



f : Respect i vely t a ken from 



Mori et al 



1 20041) . 



Webb et al. 



Kargaltsev et al. 



1 2004a 



(2001 



I 



KargaltseyfcPavlovl hooj) . 



Li et al 



Campana et al. 



1 20081) . 



DeLuca et al 



1 20051') 



DeLuca et al 



i 200a) . 



Webb et al. 



I 2004M) . 



i 200a) . 



g : The spectrum is well fitted also by a single blackbody. 



PSR Name 


(mJy) 


Ft 

(lO-i^erg/cm^s) 




77 




Cutoffc 
(GeV) 






F7 /Fx tot 


J0007+7303 


<0.006'= 


3.82±0.11 


1 


.38±0, 


.04 


4.6±0.4 


0.2 


5570±827 


4544±546 


J0030+0451 


0.6 


0.527±0.035 


1 


.22±0, 


.16 


1.8±0.4 


0.17 


454±31 


188±14 


J0205+6449 


0.04 


0.665±0.054 


2 


.09±0, 


.14 


3.5±1.4 


0.0025 


33.4±2.9 


33.4±2.9 


J0218+4232 


0.9 


0.362±0.053 


2 


.02±0, 


.23 


5.1±4.2 


0.2 


74.3t^^-j 


74.3ti3.9 


J0248+6021 


9 


0.308±0.058 


1 


.15±0, 


.49 


1.4±0.6 


0.735 


>34.2 


>34.2 


J0357+32 


<0.043= 


0.639±0.037 


1 


.29±0, 


.18 


0.9±0.2 


5.23 


1000t?^° 


1000t?°° 


J0437-4715 


140 


0.186±0.022 


1 


.74±0, 


.32 


1.3±0.7 


0.02 


18.4±2.5 


13.0±1.7 


J0534+2200 


14 


13.07±1.12 


1 


.97±0, 


.06 


5.8±1.2 


0.001 


0.295±0.026 


0.295±0.026 


J0540-6919 


0.024 


<0.833 




2 




- 


- 


<1.47 


<1.47 


J0613-0200 


1.4 


0.324±0.035 


1 


.38±0, 


.24 


2.7±1.0 


0.07 


1464tl^^^ 




J0631+1036 


0.8 


0.304±0.051 


1 


.38±0, 


.35 


3.6±1.8 


0.14 


>1350 


>1350 


J0633+0632 


<0.003'^ 


0.801±0.064 


1 


.29±0, 


.18 


2.2±0.6 


1.4 


524±179 


524±179 


J0633+1746 


<1 


33.85±0.29 


1 


.08±0, 


.02 


1.90±0.05 


0.78 


6812^136 


26871^^1 


J0659+1414 


3.7 


0.317±0.030 


2, 


.37±0, 


.42 


0.7±0.5 


0.01 


78.1113% 




J0742-2822 


15 


0.183±0.035 


1 


.76±0, 


.40 


2.0±1.4 


0.07 


>812 


>812 


J0751+1807 


3.2 


0.109±0.032 


1 


.56±0, 


.58 


3.0±4.3 


0.08 


24811^5 


24811^^ 


J0835-4510 


1100 


88.06±0.45 


1 


.57±0, 


.01 


3.2±0.06 


0.01 


1353±326 


313±75 


J1023-5746 


<0.031 


1.55±0.10 


1 


.47±0, 


.14 


1.6±0.3 


0.12 


963±173 


963±173 


J1028-5819 


0.36 


1.77±0.12 


1 


.25±0, 


.17 


1.9±0.5 


0.14 


1182±403 


1182±403 


J1044-5737 


<0.021 


1.03±0.07 


1 


.60±0, 


.12 


2.5±0.5 


0.45 


>262 


>262 


J1048-5832 


6.5 


1.73±0.11 


1 


.31±0, 


.15 


2.0±0.4 


0.08 


345lt?^f 


345lt?^f 


J1057-5226 


11 


2.72±0.08 


1, 


.06±0, 


.08 


1.3±0.1 


0.56 


1804tf^4 


111^^2 



PSRName Fr 7^ 
(mJy) (lO-i^erg/cm^s) 



J1124-5916 


0.08 


0.380±0.058 


1.43±0.33 


J1413-6205 


<0.025 


1.29±0.10 


1.32.±0.16 


J1418-6058 


<0.03^ 


2.36±0.32 


1, 


.32±0, 


.20 


J1420-6048 


0.9 


1.59±0.29 


1, 


.73±0, 


.20 


J1429-5911 


<0.022 


0.926±0.081 


1, 


.93±0, 


.14 


J1459-60 


<0.038^ 


1.06±0.10 


1, 


.83±0, 


.20 


J1509-5850 


0.15 


0.969±0.101 


1, 


.36±0, 


.23 


J1614-2230 


- 


0.274±0.042 


1. 


.34±0. 


.36 


J1617-5055 


- 


<0.3 




2 




J1709-4429 


7.3 


12.42±0.22 


1, 


.70±0, 


.03 


J1718-3825 


1.3 


0.673±0.160 


1, 


.26±0, 


.62 


J1732-31 


<0.008'= 


2.42±0.12 


1. 


.27±0. 


.12 


J1741-2054 


0.156^= 


1.28±0.07 


1. 


.39±0. 


.14 


J1744-1134 


3 


0.280±0.046 


1. 


.02±0. 


.59 


J1747-2958 


0.25 


1.31±0.14 


1, 


.11±0. 


.28 


J1809-2332 


<0.026'= 


4.13±0.13 


1, 


.52±0, 


.06 


J1811-1926 




<0.3 




2 




J1813-1246 


<0.028^ 


1.69±0.11 


1, 


.83±0, 


.12 


J1813-1749 




<0.3 




2 




J1826-1256 


<0.044^ 


3.34±0.15 


1, 


.49±0, 


.09 


J1833-1034 


0.07 


1.02±0.12 


2, 


.24±0, 


.15 


J1836+5925 


<0.0F 


e.oozbo.ii 


1, 


.35±0, 


.03 



(GeV) 



1.7±0.7 


0.01 


38.9±7.4 


34.9±6.7 


2.6±0.6 


0.43 


>287 


>287 


1.9±0.4 


0.08 


6672±3049 


6672±3049 


2.7±1.0 


0.06 


426±112 


426±112 


3.3±1.0 


0.45 


>55.0 


>55.0 


2.7±1.1 


0.52 


>269 


>269 


3.5±1.1 


0.15 


10881^^1 


10881^^1 


2.4±1.0 
- 


1.03 
- 


>958 
<4.5 


9581^3^ 
<4.5 


4.9±0.4 


0.33 


32851111 


1374+_lli 


1.3±0.6 


0.09 


240±81 


240±81 


2.2±0.3 


1.33 


>1000 


>1000 


1.2±0.2 


0.24 


977+111 

^' '-98 


977+111 

' ' -98 


0.7±0.4 


0.1 


>1030 


1030±187 


1.0±0.2 


0.02 


26.9^4^3^ 


26.9^4^3^ 


2.9±0.3 


0.33 


295lt^^^ 
<11.25 


1316+2^2 

10iO_226 

<11.25 


2.9±0.8 


0.20 


175±59 
<11.25 


175±59 
<11.25 


2.4±0.3 


20.7 


2834±1398 


2834±1398 


7.7±4.8 


0.01 


15.3±1.9 


15.3±1.9 


2.3±0.1 


2 


13065+"'*'^'* 

iOUDO_4958 


105201^231 



PSR Name 




F, 






77 




Cutoffc 


Eff^ 


F^/Fxnt 


F^/Fjftot 






(mJy) 




cm?s) 








(GeV) 








J1846+0919 


<0.004 


0.358±0 


035 


1 


60±0 


19 


4.1±1.5 


2.1 


>123 


>123 


J1907+06 


0.0034'' 


2.75±0 


13 


1 


84±0 


08 


4.6±1.0 


0.30 


700±260 


700±260 


J19524 


-3252 


1 


1.34±0 


07 


1 


75±0 


10 


4.5±1.2 


0.02 


38.3±5.3 


35.2±3.4 


J19544 


-2836 


<0.004 


0.975±0 


068 


1 


55±0 


14 


2.9±0.7 


0.39 


>267 


>267 


J19574 


-5036 


<0.025 


0.227±0 


020 


1 


12±0 


28 


0.9±0.2 


5.6 


>76.2 


>76.2 


J19584 


-2841 


<0.005^ 


0.846±0 


069 





77±0 


26 


1.2±0.2 


1.6 


539±187 


539±187 


J20214 


-3651 


0.1 


4.70±0 


15 


1 


65±0 


06 


2.6±0.3 


0.07 


01 oq+343 


' °"J-449 


J20214 


-4026 


<o.oir 


9.77±0 


18 


1 


79±0 


03 


3.0±0.2 


2.2 


22061±7381 


22061±7381 


J20324 


-4127 


0.05^= 


1.12±0 


12 





68±0 


38 


2.1±0.6 


0.64 


2636±790 


2636±790 


J20434 


-2740 


7 


0.155±0 


027 


1 


07±0 


55 


0.8±0.3 


0.09 


747+217 

' ^ ' —409 


747+217 

' ^ ' —409 


J2055 


-f25 


<0.106 


1.15±0 


07 





71±0 


19 


1.0±0.2 


5.4 




3010^^^^^ 


J2124-3358 


1.6 


0.276±0 


035 


1 


05±0 


28 


2.7±1.0 


0.05 


410+107 


900+74 


J22294 


-6114 


0.25 


2.20±0 


08 


1 


74±0 


07 


3.0±0.5 


0.025 


42.91^:? 




J2238+59 


<0.007'= 


0.545±0 


059 


1 


00±0 


36 


1.0±0.3 


0.52 


>121 


>121 



Table 3:: 7-ray spe c tra o f the pulsars. A broken powerlaw spectral shape is assumed for all the pulsars and the val ues are taken fro m lAbdo et al.l 1 2010al ) 



Saz Parkinson et al 



1 2OI0I ). The gamma-ray flux is above 100 GeV. T he 4 sources with an upper limit flux are taken from 



radio flux densities (at 1400MHz) are taken from 



Abdo et al 



i201oi: 


Saz Parkins 


on et al. 


1. See 


Watters et al. 


1 


2009 


)• 



Bird et al 



I 20091) (see section EAJI. The 



(2010). All the errors are at a 90% confidence level. 



b : taken from 



Abdo et al. 



c : taken from I Ray et al 



I gOlOel 'l 



1 20101 ) 



|2009|). Here the errors arc not reported. 
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Fig. 1. — : E-hx diagram for all pulsars classified as type 2 and with a clear distance estimation, assuming fx = l (see Equation 
5). Black: radio-quiet pulsars; red: radio-loud pulsars; blue: millisecond pulsars. The linear best fit of the logs of the two 
quantities is shown. 
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Fig. 2. — : E-hj diagram for all pulsars classified as type 2 and with a clear distance estimation, assuming f^ = l (see Equation 
5). Black: radio-quiet pulsars; red: radio-loud pulsars; blue: millisecond pulsars. The double linear best fit of the logs of the 
two quantities is shown. 
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Fig. 3. — : log(FT,/Fx) histogram. The step is 0.5; the radio-loud (and millisecond) pulsars are indicated in grey and the 
radio-quiet ones in black. Only high confidence pulsars (type 2) have been used for a total of 24 radio-loud and 9 radio-quiet 
pulsars. 
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Fig. 4. — : E-F^,/Fx diagram. Green: IBIS pulsars; black: radio-quiet pulsars; red: radio-loud pulsars; blue: millisecond 
pulsars. The triangles are upper and lower limits, the squares indicate pulsars with a type 1 X-ray spectrum (see Table [2]l and 
the stars pulsars with a high quality X-ray spectrum The dotted line is the combination of the best fitting functions obtained 
for Figure [T] and [2] with the geometrical correction factor set to 1 for both the X and 7-ray bands. 
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Fig. 5. — : E-Fj/Fx diagram for high confidence pulsars only (type 2). Green: IBIS pulsars; black: radio-quiet pulsars; red: 
radio-loud pulsars; blue: millisecond pulsars. The triangles are upper limits. 
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Fig. 6. — : Left: Characteristic Age-F^ /Fx diagram. Right: SNR Age-F-y/Fx diagram. Green: IBIS pulsars; black: radio-quiet 
pulsars; red: radio-loud pulsars. Triangles are upper limits, squares are pulsars with a type 1 X-ray spectrum while stars are 
pulsars with a type 2 X-ray spectrum. 



